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Abstract. We derive the maximum accelerated particle en- 
ergy from first-order and second-order Fermi acceleration 
at nonrelativistic and relativistic shocks for explosions tak- 
ing place in a uniform surrounding medium. Second-order 
stochastic processes in relativistic flows are shown to be ca- 
pable of accelerating cosmic rays to ultra-high energies. Cos- 
mic rays above the knee of the cosmic ray spectrum can be 
accelerated by the second-order Fermi mechanism in rela- 
tivistic flows, such as those occuring in gamma-ray bursts 
and unusual supernovae like SN 1998bw. 



1 Introduction 

Although there is a consensus that cosmic rays with ener- 
gies less than the knee energy at w 3 x 10 15 eV are accel- 
erated by supernova remnant shocks, the origin of higher en- 
ergy cosmic rays is unknown. Proposed sources of cosmic 
rays above the knee include galactic sources such as micro- 
quasars and pulsars, active galactic nuclei (Protheroe and Sz- 
abo, 1992), merger shocks in galaxy clusters (Kang Ryu, and 
Jones, 1996), dormant AGNs (Levinson and Boldt, 2001), 
and gamma-ray bursts (GRBs) (Milgrom and Usov, 1996; 
Dermer and Humi, 2001). Proposed sources of ultra-high 
energy ( ^ 10 19 eV) cosmic rays (UHECRs) include the ter- 
mination shocks at the lobes of FRII jets (Rachen and Bier- 
mann, 1993), GRBs (Waxman, 1995; Vietri, 1995; Dermer, 

2000) , and exotic particles, such as monopoles and Z-bursts 
(for reviews, see Nagano and Watson (2000); Weiler (2001)). 

Here we consider particle acceleration at nonrelativistic 
and relativistic shocks. We show that the sources of cos- 
mic rays near and above the knee of the cosmic-ray spec- 
trum could result from particle acceleration at the relativis- 
tic shocks produced by the baryon-dilute outflows of fire- 
ball transients and GRBs (Dermer, 2000; Dermer and Humi, 

2001) . The maximum energy that can be reached by a par- 
ticle accelerated at a shock due to constraints upon diffu- 
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sive transport and available time is derived, extending ear- 
lier treatments for nonrelativistic shocks by, e.g., Lagage and 
Cesarsky (1983), Drury (1983), and Baring et al. (1999). 

In §2, the maximum particle energy is derived for the ex- 
tremes of quasi-parallel (q-||) and quasi-perpendicular (q-_L) 
shocks (Jokipii, 1987; Baring et al., 1999). We find that the 
maximum particle energy is not greatly increased in the latter 
case due to a geometric constraint that applies in cases where 
a large scale toroidal magnetic field is absent. In §3, we con- 
sider the highest energy that particles can be accelerated by 
shock Fermi processes at a relativistic external shock, follow- 
ing Gallant and Achterberg (1999). In §4 we show that the 
highest energies that can be reached by Fermi processes in 
supernova or GRB-type explosions result from second-order 
gyroresonant stochastic acceleration in the shocked fluid of a 
relativistic blast wave (Waxman, 1995; Rachen and Meszaros, 
1998; Schlickeiser and Dermer, 2000; Dermer and Humi, 
2001). This maximum particle energy E can exceed w 10 20 
e V in GRB blast waves for reasonable parameter values. Cos- 
mic rays with energies above the knee of the cosmic ray spec- 
trum is proposed to originate from second-order Fermi accel- 
eration at relativistic shocks (§5). 



2 Nonrelativistic Shocks 

We treat first-order Fermi acceleration of relativistic nonther- 
mal particles with v' w c. Our notation is that j3(x)c = (3c 
is the flow speed, and the blast-wave evolution is described 
by momentum P = (3T. Primes refer to quantities in the 
comoving shock frame. In the primed frame stationary with 
respect to the shock, the upstream (— ) flow approaches with 
speed u_ = /3_c = 4/3c/3 and the downstream (+) flow re- 
cedes with speed u + = (3 + c = /3c/3. The quantity u = 
u- — u+ — (3c is the speed of the shocked fluid and u_ is 
the speed of the shock as measured in the stationary upstream 
frame. 

Let p' — \/~f' 2 — 1 represent the dimensionless momen- 
tum of a particle with Lorentz factor 7' in the comoving 
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shock frame. Then p' FI = Ap' /t cyc , where the cycle time 
t C yc is given by the diffusive properties of the upstream and 
downstream regions. The average change in particle momen- 
tum over a complete cycle for relativistic test particles with 
p' > 1 is Ap' = 40p'/3, when < 1 (Schlickeiser, 1984; 
Gaisser, 1990), provided that there are scatterers in the up- 
stream and downstream flow to isotropize the particles. In a 
one-dimensional flow, 
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where x = u-/u+ is the compression ratio and the diffu- 
sion coefficient n± — X±v'/3 = rj±r^v'/3 = r]±r^ ± p'v' /3, 
where r^± = r^ ± p' = mc 2 p' /qB± is the Larmor radius for 
a particle of mass m and charge q. Here we write diffusion 
coefficients in terms of the parameters n± that give the parti- 
cle mean-free-paths scaled to the values implied by the Bohm 
diffusion limit evaluated for the local magnetic field B. 
In q-|| shocks, B-/B+ = 1. Hence 
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The last relation defines the Bohm-diffusion-limit maximum 
acceleration rate for q-|| shocks, assuming that the Bohm dif- 
fusion limit sets a lower limit to the diffusivity of a medium. 

In q-_L shocks, a scattering event shifts a particle in the 
x direction by a mean distance t"l; therefore A — ► n,. The 
effective particle drift speed is reduced by a factor rf 1 be- 
cause particles are confined within the gyro-radius size scale 
n,. Hence k± = r^v' /3n or k± = K||/(l + r/ 2 ) (Jokipii, 
1987), and 
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where 77 = min(ri_ ,??+)• 

According to Jokipii (1987), to maintain an isotropic par- 
ticle distribution requires rj v'/u ~ 1//3 for relativistic 
particles. If shock drift takes particles to regions of smaller 
obliquenesses, then the acceleration rate is reduced (Jokipii, 
1987; Baring et al., 1999). Particles will drift along a q- 
_L shock to scatter in a region of smaller obliqueness when 
ryrL — /x £ #1/2 , where &1/2 is a characteristic angle over 
which the shock obliquity changes. If the surrounding mag- 
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netic field is homogeneous, ry_ 
condition 77 < 
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y/x/tx- implies that r\ 0~ 
The energy gain rates are more conveniently written in 
terms of the differential distance dx = s T s cdt' — > s cdt' = 
U-dt' = u_dt that the shock advances during the differential 
time element dt, as measured in the stationary frame of the 
external medium. From the above results, we find that the 
maximum energy gains per unit distance for quasi-parallel 



and quasi-perpendicular nonrelativistic shocks are, respec- 
tively, 

^-\ lmax ~0 qB - ■ ^\^ max ^0 2 /\ B . . (5) 

The maximum particle energy in nonrelativistic supernova 
shocks is derived by considering a uniform medium with 
density no that carries a homogeneous external magnetic field 
with intensity B- = B n = 10^ 6 B^a G. Shock-wave evo- 
lution is assumed to follow the behavior of the coasting and 
adiabatic Sedov solutions. The evolution of the blast wave 
momentum P{x) = 0{x)T{x) in the adiabatic (Sedov) regime 
is described by the function 
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(Dermer and Humi, 2001), where P = ^Tg - 1 = /? r is 
the initial blast-wave Lorentz factor that defines the baryon 
loading. The deceleration radius 
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(Meszaros and Rees, 1993) is defined in terms of the di- 
rectional energy release m of the explosion, expressed in 
units of Solar rest-mass energy per 4ir sr. Thus m Q = 1 
dE /dfl = Mqc 2 /4tt = 1.4 x 10 53 ergs sr" 1 . For a non- 
relativistic explosion, x d = Is, the Sedov length (e.g., Sari 
et al., 1996), where £ s = T 2 Q /3 x d . 

The maximum particle energy for a particle accelerated 
from rest starting near the location of the explosion is found 
by integrating the energy gain rate per unit distance over the 
evolution of the blast wave radius. For q-|| and q-_L nonrela- 
tivistic shocks, we obtain 
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respectively, where Ze is the particle charge. 



3 Relativistic Shocks 

It is simplest to work in the frame of the upstream medium 
to calculate the maximum energy of particles accelerated by 
a relativistic external shock. In the first cycle, a particle in- 
creases its energy by a factor T 2 but, due to kinematics of 
escape and capture, the particle increases its energy by only 
a factor w 2 in successive cycles, as shown by Gallant and 
Achterberg (1999). This is because a particle is captured by 
the advancing relativistic shock when the particle is deflected 
by an angle 9d <~ 1/T S , so that the upstream cycle time is 
t u ~ (r s w^_) _1 . Here the particle gyration frequency in 
the upstream medium is lj±- = v/tl-, where v is the parti- 
cle speed and the particle deflection is determined by £>j__, 
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the transverse component of the magnetic field in the up- 
stream region. As a result, the nonthermal particle distribu- 
tion function at a relativistic shock is highly anisotropic (Kirk 
and Schneider, 1987). As shown by Gallant and Achterberg 
(1999), the cycle time is dominated by the upstream transit 
time, 

Because the energy increases by a factor = 2 during a sin- 
gle cycle, the momentum gain rate following the first cycle 
in relativistic shock acceleration is 
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provided T s » 1. If the particles are captured in the first cy- 
cle with Lorentz factor 7, then the energy from the first cycle 
of Fermi acceleration in the stationary frame is w T 2 jmc 2 
(Vietri, 1995). After integrating the energy gain rate, we find 
that the maximum energy that can be achieved by relativistic 
first-order shock acceleration at an external shock is therefore 
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where Am p is the particle mass, and T = 300r 3 oo- From 
equation (1 1), we see that it is difficult to accelerate particles 
to energies larger than a factor <~ 2 1 ^ 3 T^ 3 over Ej_^ max , 
given by equation (9), unless a pre-existing energetic particle 
distribution is found in the vicinity of the explosion. 



4 Gyroresonant Stochastic Acceleration in Shocks 

We derive the particle acceleration rate through stochastic 
gyroresonant processes from the expression p' F11 = Ap'/ti SO , 
where ti SO is the pitch-angle isotropization time scale in the 
comoving fluid frame. For relativistic hard sphere scatter- 
ers, the fractional change in momentum over time period 
t iso is Ap'/p' = 4p^/3 (Gaisser, 1990). The quantity pac, 
which represents the dimensionless momentum of the scat- 
tering centers, reduces to the Alfven speed in the weakly tur- 
bulent quasilinear regime. For relativistic shocks, the maxi- 
mum value that p\ can take is given by p 2 A = 2e B /3, as can 
be shown by considering the relativistic shock jump condi- 
tions (Blandford and McKee, 1976). The downstream mag- 
netic field B + = max(xr J B_ L _, J B,Vr 2 -T ), where B* = 
{8irxn m p c 2 e B ) 1/2 = 0.39(%n o e B ) 1/2 G. The first term 
represents the compression of the upstream transverse mag- 
netic field, and the second term defines the downstream field 
in terms of the equipartition field B eq (G) — 0.39(e B no) 1 ^ 2 
Vr 2 — T through the parameter e B . For this convention, 
X = 4 for strong shocks whether or not the shocks are rela- 
tivistic. 

The particle pitch angle changes by SB /B during one gy- 
roperiod t gyr = r^+jc, where B now refers to the mean 
downstream magnetic field in the shocked fluid. Thus ti SO = 
t gyr /(SB/B) 2 , noting that particles diffuse in pitch angle. 



Pitch angle changes due to gyroresonant interactions of par- 
ticles with resonant plasma waves are described by the rela- 
tion (SB/B) 2 w kW(k)/U B , where U B = B 2 /8n, and the 
turbulence spectrum W(k) = Wo(k/k min )~ v for k min < 
k < k max . The index v = 5/3 for a Kolmogorov turbu- 
lence spectrum, and v = 3/2 if the Kraichnan phenomenol- 
ogy is used. The resonant wavenumber is normally assigned 
through the resonance condition oj — k^v^ = ifl/j, but here 
we proceed by employing the simple resonance assumption 
k — > 1/tl+ (Biermann and Strittmatter, 1987). Assuming 
isotropy of forward and backward-moving waves gives the 
normalization Wo = £Ub (v — 1)/ 2fc TO ,„, where £ is the ratio 
of plasma turbulence to magnetic field energy density. Hence 
tul - c ^ v - l)(rL+k mm r- 1 /(2r L+ ), and 
., 2 



l\v-l 



(12) 



Pfii - ~Pa £{v ~ 1) (-3-) {rl+kminP' 

(e.g.,Dermeret al. (1996)). 

The term r^ + k min p' in the parentheses of equation (12) 
gives the comoving gyroradius in units of the inverse of the 
smallest turbulence wavenumber k min found in the shocked 
fluid. If k min ~ 1/ A', where A' is the blast wave width, 
then the relation r£ + fc m i n p' ^ 1 is an expression of the con- 
dition that the particle gyroradius is less than the size scale 
of the system (Hillas, 1984). When this condition is satisfied, 
PFii, m «»l^c/r£ + S2eB C /(3r£ + ). 

By integrating the energy gain-rate for an adiabatic blast 
wave that evolves according to equation (6), we obtain the 
maximum energy-gain rates due to stochastic Fermi acceler- 
ation for nonrelativistic (T — 1 <C 1) and relativistic (T >• 1) 
shocks, given by 

2 1 / 2 E' 



dE' , 
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and 
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(14) 



respectively. Here /a sets the scale for the smallest turbu- 
lence wave numbers in comparison with the inverse of the 
comoving size scale x/T of the blast wave width. From hy- 
drodynamic considerations, /a = 1/12. 

For an adiabatic blast wave that explodes in a uniform sur- 
rounding medium, the nonrelativistic expression (13) can be 
integrated to give the maximum energy 



,FII,NR — Kjsi qB^fAXdPo 



(15) 



where 
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/a 
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Here /3 = 10 2 (3-2, and I 2 n is an integral that takes a max- 
imum value ~ 6. 
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The maximum energy of particles accelerated by second- 
order Fermi acceleration in a nonrelativistic blast wave is de- 
fined by the quantity 

qB*f A x d p £* 8 x 10 18 Z e ^ /2 (n ) 1/6 /A/3-2m^ /3 eV , (17) 

though the actual maximum energy is much less than this 
value when /3 <C 1, as can be seen from equation (16). 
The relativistic expression (14) can be integrated to give 

E max , FII , E R = e B&V-l) ^. qB j AXdTo (lg) 

The maximum energy of particles accelerated by second- 
order Fermi acceleration in a relativistic blast wave occurs 
near the deceleration radius x w xj. The basic scaling for 
this maximum energy is given by the term 

qB*f A x d T = 7.7xl0 20 Z e ^ /2 ny 6 / A (m Q r ) 1/3 cV . (19) 

Note the sensitive dependence of the maximum energy given 
byequation(18)onthefactores^//A- If es£//A > 1 when 
the blast wave is near the deceleration radius, then stochastic 
acceleration processes in relativistic blast waves can in prin- 
ciple accelerate particles to ultra-high energies. 

5 Discussion 

The maximum particle energy that can be achieved through 
first-order Fermi acceleration of a particle initially at rest, for 
an adiabatic blast wave in a uniform surrounding medium, is 
given by 

E maxA c i^ZB^Ti 7 ^)^ eV . (20) 

Here we use the q-_L result for nonrelativistic blast waves, 
which gives the largest energies. The maximum particle en- 
ergy that can be achieved through second-order Fermi accel- 
eration in a relativistic blast wave is given by 

E max , 2 * 8 x l0 20 K v Ze]i 2 ny 6 f A {mQT f^ eV , (21) 

where K v = [2 3 / 2 e B £/9fA] 1/( - 2 ~ v) - The strong dependence 
of the maximum energy on (3q makes second-order processes 
unimportant in nonrelativistic flows when /3rj <C 0.1. As the 
flow speeds become marginally relativistic, first-order pro- 
cesses may accelerate particles to sufficiently high energies 
that second-order acceleration then starts to become impor- 
tant (Schlickeiser, 1984). At relativistic shocks, the second- 
order process can easily dominate. 

Spectral models of GRBs within the external shock sce- 
nario (Meszaros and Rees, 1993; Chiang and Dermer, 1999) 
imply r 30 o ~ 1 and predict the existence of clean and dirty 
fireballs which have not yet been detected (Dermer et al., 
1999). Mildly relativistic outflows are deduced from radio 
observations of the Type Ic SN 1998bw associated with GRB 
980425 (Weiler et al., 2000; Kulkarni et al., 1998), so that 
hadronic cosmic ray acceleration to energies above the knee 
of the cosmic ray spectrum is possible from Type Ic SNe. 



First-order Fermi processes in relativistic flows can accel- 
erate particles to energies above the knee of the cosmic -ray 
spectrum. The smooth cosmic ray spectrum between the 
knee and the ankle suggests that the cosmic rays with en- 
ergies above the knee originate from sources that accelerate 
particles to energies £ 10 19 eV. The second-order process 
in relativistic shocked fluids can accelerate particles to ultra- 
high energies to form both the galactic halo cosmic rays with 
energies between the knee and the ankle, as well as the meta- 
galactic UHECRs. The GeV cosmic rays detected near Earth 
are probably a mixture of accelerated particles produced by 
many sources, but the cosmic rays formed near and above 
the knee may be produced by only one or a few energetic ex- 
plosions (Erlykin and Wolfendale, 2000). GRBs are a prime 
candidate for the origin of these cosmic rays (Milgrom and 
Usov, 1996; Dermer and Humi, 2001). 
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